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O ' Abstract. We present the asymptotic formula for the Wigner 9j-symbol, valid when 

0^ . all quantum numbers are large, in the classically allowed region. As in the Ponzano- 
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Rules are presented for converting spin networks into the associated geometrical figures. 
The amplitude is expressed as the determinant of a 2 x 2 matrix of Poisson brackets. 
The 9j-symbol possesses caustics associated with the fold and elliptic and hyperbolic 
umbilic catastrophes. The asymptotic formula obeys the exact symmetries of the 9j- 
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^ ! 1. Introduction 



The asymptotic behavior of spin networks has played a significant role in simplicial 



approaches to quantum gravity. Indeed, the field began with the observation that the 
Ponzano-Regge action (1968) for the semiclassical 6j-symbol is identical to the Einstein- 



Hilbert action of a tetrahedron in 3-dimensional gravity in the Regge formulation (Regge, 
1961; see also Williams and Tuckey 1992 and Regge and Williams 2000). More recently, 
semiclassical expansions have been used to study the low-energy or classical limit of 
quantum gravity as well as to derive quantum corrections to the classical theory. 
Asymptotic studies in this area have included treatments of the lOj-symbol (Barrett 
and Williams 1999, Baez et al 2002, Barrett and Steele 2003, Freidel and Louapre 
2003), amplitudes in the Freidel-Krasnov model (Conrady and Freidel 2008), LQG 
fusion coefficients (Alesci et al 2008), and the EPRL amplitude (Barrett et al 2009). 
In addition, the venerable 6j-symbol and Ponzano Regge (1968) formula continue to 
receive attention (Roberts 1999, Barrett and Steele 2003, Freidel and Louapre 2003, 
Gurau 2008, Charles 2008, Littlejohn and Yu 2009, Depuis and Livine 2009, Ragni et al 
2010), not to mention the g-deformed 6j-symbol (Nomura 1989; Taylor and Woodward 
2004, 2005). 

In this article we present the generalization of the Ponzano-Regge formula to the 
Wigner 9j-symbol, as well as some material relevant for the asymptotics of arbitrary spin 
networks. The Ponzano-Regge formula (Ponzano and Regge 1968) gives the asymptotic 
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expression for the Wigner 6j-symbol when all quantum numbers are large. The 9j- 
symbol is the next most complicated spin network after the 6j-symbol, with features 
that are found in all higher spin networks. In this article we present only the asymptotic 
formula itself for the 9j-symbol and some salient facts surrounding it. We defer a 
derivation and deeper discussion of the formula to a subsequent publication. 

Our derivation has quite a few steps, and some of them at this point are supported 
by numerical evidence only. Thus, we do not now have a rigorous derivation of our 
result. We believe it is correct, however, on the basis of direct numerical comparisons 
with the exact 9j-symbol, the fact that our formula obeys all the symmetries of the exact 
9j-symbol, and the plausibility and numerical support for the conjectures involved in 
the parts of the derivation currently lacking proofs. The proofs do not seem difficult, 
and we hope to fill in the gaps in our future work. 

Although most of the papers cited above have dealt with the asymptotics of specific 
spin networks, usually there are special values of the angular momenta that are used. For 
example, the lOj-symbol involves balanced representations of 5*0(4), which means that 
some pairs of j's are equal, while the 9j-symbols that appear in LQG fusion coefficients 
have two columns in which one quantum number is the sum of the other two. In addition, 
j's are sometimes set equal because this is regarded as the most interesting regime from 
a physical standpoint. 

As a result, the spin networks that have been studied tend to fall on caustics where 
the asymptotic behavior is not generic. At such points, the value of the spin network 
(the wave function) is not oscillatory in a simple sense, instead it has the form of a 
diffraction catastrophe (Berry 1976). In addition, the wave function scales as a higher 
(less negative) power of the scaling parameter (effectively, 1/H). This type of behavior 
has been noted in several places in the quantum gravity literature, although as far as 
we can tell no one has noted that it is related to standard caustic and catastrophe 
types. In this article we give a rather complete picture of the 9j-symbol for all possible 
parameters in the classically allowed region, including all phases and Maslov indices. 
We also indicate the subsets upon which the behavior is nongeneric and described by 
various types of caustics. We believe that this is the first time that such information 
has been available for any spin network more complicated than the 6j'-symbol. 

Another reason for interest in the 9j-symbol is that it is the nontrivial part of the 
Clebsch-Gordan coefficient for SO (4). 

Basic references on the Wigner 9j-symbol include Edmonds (1960), Biedenharn 
and Louck (1981ab) and Varshalovich et al (1981). Recent work on the 9j-symbol has 
included new asymptotic forms when some quantum numbers are large and others small 
(Anderson et al 2008, 2009). We also note the use of SU{2) spin networks in quantum 
computing (Marzuoli and Rasetti 2005). 

In Sec.[2]we present the asymptotic formula for the 9j-symbol and draw comparisons 
with the Ponzano-Regge formula to introduce its geometrical content. A detailed 
explanation of the notation follows in later sections. In Sec. [3]we present general rules for 
converting spin networks into surfaces composed of oriented edges and oriented triangles, 
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and illustrate them for the 9j-symbol. In Sec. S]we explain how the geometrical objects 
(pieces of oriented surfaces) corresponding to the 9j-symbol can be constructed in 3- 
dimensional space. In Sec.[5]we explain the configuration space of the 9j'-symbol and the 
classically allowed subset thereof. In Sec. [6] we define the amplitude of the asymptotic 
formula and discuss the manifolds (the caustics) upon which it diverges as well as the 
diffraction catastrophes that replace the simple asymptotic form in the neighborhood 
of the caustics. In Sec. [7] we explain the phase of the semiclassical approximation, a 
generalization of the Ponzano-Regge action that requires careful definitions of dihedral 
angles. In Sec. [8] we show that the asymptotic formula correctly obeys the symmetries 
of the 9j-symbol. Finally, in Sec. [9] we present some comments and conclusions. 

2. The Asymptotic Formula 

The asymptotic expression for the 9j-symbol is 



where A\p are positive amplitudes, S\$, are phases, and each term is roughly similar to 
the single term in the Ponzano-Regge formula for the 6j-symbol. The right hand side is 
the leading term in an asymptotic expansion in powers of 1/k of the 9j-symbol when all 
nine j's are scaled by a positive factor k that is allowed to go to infinity (k plays the role 
of 1/h in the asymptotic expansion). The fc's are suppressed in ([1]), but the expression 
on the right scales as 1/k 3 . Equation (CQ) applies only in the classically allowed region. 
We do not present the analog of ([1]) in the classically forbidden region. 

Equation (OQ) breaks down near caustics, where the 9j-symbol scales with a higher 
(less negative) power of k than 1/k 3 . In the neighborhood of caustics, the 9j-symbol is 
approximated by diffraction catastrophes, including the fold and hyperbolic and elliptic 
umbilic. These are discussed more fully in Sec. El 

To explain the meaning of ([I]) some analogies with the Ponzano-Regge formula 
for the 6j-symbol are useful. In the classically allowed region, the Ponzano-Regge 
formula associates a given 6j-symbol with a real tetrahedron whose six edge lengths are 
J% — ji + 1/2, where the six j's are those appearing in the 6j-symbol. More precisely, 
there are two tetrahedra, related by spatial inversion, that is, time-reversal. Except for 
flat configurations, the two tetrahedra are not related by proper rotations in SO (3). We 
recall that time-reversal, not parity, inverts the direction of angular momentum vectors. 
The two tetrahedra correspond to the two stationary phase points of the 6j-symbol, 
which make contributions to the asymptotic expression that are complex conjugates of 
each other. The result is the real cosine term in the Ponzano-Regge formula. One can say 
that semiclassically the 6j-symbol is a superposition of two amplitudes, corresponding 
to a tetrahedral geometry and its time-reversed image, that produce oscillations in the 
result. 




(1) 
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We shall use lower case j's for quantum numbers, and capital J's for the lengths of 
the corresponding classical vectors. These are always related by Jj = jj + 1/2. The 1/2 
is a Maslov index (Maslov and Fedoriuk 1981, Mischenko et al 1990, de Gosson 1997), 
and the manner in which it arises in this context is explained in Aquilanti et al (2007). 

In the case of the 9j-symbol in the classically allowed region, there are four 
geometrical figures associated with a given set of nine j's, consisting of two pairs related 
by time-reversal. The four geometrical figures correspond to the four real stationary 
phase points of the 9j-symbol. Each pair of figures is associated with an "admissible" 
root (defined momentarily) of a certain quartic equation. There are two admissible roots 
in the classically allowed region, labeled 1 and 2, corresponding to the two terms in ([TJ). 
Each trigonometric term in (JI]) consists of an exponential and its complex conjugate, 
corrresponding to a geometrical figure and its time-reversed image. One can say that 
semiclassically the 9j-symbol is a superposition of four amplitudes corresponding to four 
geometries, consisting of two pairs of a geometry and its time-reversed image. We now 
explain these geometries and how they are specified by the nine j's that appear in the 
symbol. 

3. Triangles, Orientations and Geometries 

The 9j-symbol specifies the lengths Jj = jj + 1/2 of nine classical angular momentum 
vectors Jj but not their directions. Therefore we inquire as to how the directions may 
be determined, and geometrical figures constructed out of the resulting vectors. 

Actually, it is convenient to double this set and speak of 18 classical vectors Ji, 
3'i, % — 1, . . . , 9. A doubling of this kind was introduced by Roberts (1999), who gave 
a highly symmetrical way of writing the 6j-symbol as a scalar product in a certain 
Hilbert space. Although Roberts only worked with the 6j-symbol, his method is easily 
generalized to an arbitary spin network. Ponzano and Regge (1968) also gave hints that 
doubling of angular momentum vectors are important in the asymptotic analysis of spin 
networks. 

We now describe rules that take an arbitary spin network (with at most trivalent 
vertices) and transcribe it into relations among a doubled set of classical angular 
momentum vectors, defining a set of oriented triangles and oriented edges of a 
geometrical figure. We exemplify these rules only in the case of the 9j-symbol, but 
they are easily applied to any spin network. The reader may find it illuminating to 
apply our rules to the 6j-symbol, starting with the usual spin network (the Mercedes 
graph). Figure [1] illustrates the spin network of the 9j-symbol. See also Fig. 18.1 of 
Yutsis et al (1962). 

Each edge of the spin network, labeled by jj, is associated with two classical angular 
momentum vectors Jj and 3\ that are required to satisfy 

ijji = ij:i = ^=j- + i/2 (2) 
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and 



j 4 + j; = o. 



(3) 



Vectors Jj and have the same length and point in opposite directions. 

Each vertex of the spin network, where three edges meet, corresponds to three 
vectors that add to zero. The three vectors are associated with the three edges. If the 
arrow on an edge ending at the vertex is pointing away from the vertex, then the angular 
momentum vector is unprimed; if it is pointing toward the vertex, then the vector is 
primed. This rule applied to Fig. [I] gives 



These are a set of classical triangle relations, one for each vertex of the spin network. 
In the case of the 9j-symbol, they are obviously related to the rows and columns of the 
symbol. 

Although the vector addition in PJ) is commutative, we agree to write the vectors 
in each equation in counterclockwise order (around the vertex of the spin network) for a 
vertex with + orientation, and in clockwise order for a vertex with — orientation, modulo 
cyclic permutations. Thus the ordering of the vectors is the same as the ordering of the 
columns of the 3j-symbol implied by the vertex of the network. 

This ordering is used to define a set of oriented triangles. We take the three vectors 
of any one of the equations (JI]) and place the base of one vector at the tip of the preceding 
one, to create the three edges of a triangle. In this process we parallel translate the 
vectors (in R 3 ) but do not rotate them. The triangle is given an orientation (a definition 
of a normal) by taking the cross product of any two successive vectors defining the edges. 
For example, the normal to the 123-triangle is Ji x J 2 , and that of the l'4'7'-triangle is 
J' x x J' 4 , which, in view of d2J), is the same as Ji x J 4 . 

Next, we take the triangles and displace them so that the edge Jj of one triangle is 
adjacent to the edge J[ of another triangle. In this process, the triangles are displaced 
but not rotated. If we do this with the six triangles defined by (HI) in the case of the 
9j-symbol, we find that six pairs of edges can be made adjacent, as illustrated by the 
central six triangles of Fig. [2 In this "central region" six pairs of vectors J, and J' { are 
adjacent for i = 1,2, 5, 6, 7, 9. There is some arbitrariness in choosing which six pairs 
of edges will be made adjacent. If we wish that the remaining edges i = 3, 4, 8 also be 
paired, we can duplicate three of the triangles and attach them to the periphery of the 
central region, as illustrated in Fig. [2j This amounts to a kind of "analytic continuation" 
of the central region. 

Figure [2] is highly schematic. In general, the triangles are not equilateral, the 
surface that is formed by attaching them together is not planar, and the triangles may 
fold under one another. 

The central region in Fig. [2] is a piece of an oriented surface, that is, all the normal 
vectors (by our convention) are pointing on the same side. In the case of the 6j-symbol, 



Jx + J 2 + J 3 = 0, 

J 4 + J 5 + Je = 0, 
J 7 + J 8 + J 9 = 0, 




o. 

0. 



0. 



(4) 
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our rules produce a closed surface (the usual tetrahedron), with normals all pointing 
either outward or inward (time-reversal converts one into the other). In the case of 
the 9j-symbol, the surface is not closed. There is some suggestion that this surface 
represents a triangulation of RP 2 but for this article we shall view it as living in R 3 . 

Finally, we orient each edge by choosing the direction of the vector Jj (not JQ. 

We will be interested in finding solutions {Jj, J^, % — 1, . . . , 9} of (j2J), © and (j4]), 
modulo overall proper rotations (in 5*0(3)). That is, although we do not rotate vectors 
or faces when forming our surface with oriented faces and egdes, we are allowed to rotate 
the whole surface once completed. 

We notice that if {Jj,J^,z = 1,...,9} is a solution of these equations, then the 
time-reversed set {— J i; — J-, i = 1, . . . , 9} is also a solution. If we apply our rules for 
converting vectors into a surface, we will find in general that the time-reversed set 
produces a different surface (not equivalent under SO (3)). We apply time- reversal only 
to the vectors, not the rules; for example, the ordering of the time-reversed vectors is 
the same as the original vectors. The central six triangles of the time-reversed surface 
are illustrated in Fig. [3j 

To visualize the surfaces in Figs. [2] and [3J we may imagine that the central region 
of Fig. [2] bulges out of the paper, like the northern hemisphere of a sphere (whether 
it does or not depends on the parameters, but this is one possibility). Then the time- 
reversed surface in Fig. [3] bulges into the paper, since spatial inversion is equivalent, 
modulo S0(3), to reflection in a plane. Then the central region of Fig. [2] can be glued 
to the time-reversed surface in Fig. El bringing edge J' 3 adjacent to edge — J3, etc, and 
producing a surface homeomorphic to S 2 . This is the hexagonal bipyramid constructed 
by Ponzano and Regge (1968). The conventional normals are pointing outward in the 
northern hemisphere, and inward on the southern. As noted by Ponzano and Regge, 
this bipyramid is bisected by three planes passing through a common line, namely the 
"axis" of the sphere, which cut the bipyramid into three pairs of congruent tetrahedra. 
These correspond to the three 6j-symbols in the representation of the 9j-symbol as a 
sum over products of 6j-symbols (see Edmonds (1960) Eq. 6.4.3), in which the variable 
of summation is the common edge of the tetrahedra (the axis of the sphere). 

4. Finding the Vectors 

To find a solution of (j2J), © and P| we notice that all 18 vectors are determined 
if only four of them, {Ji, J2, J4, J5} are given. We let G be the 4x4 Gram matrix 
constructed out of these vectors, that is, the 4x4, real symmetric matrix of dot products 
of these vectors among themselves. Of the ten independent dot products, eight can be 
determined from the given lengths Jj, i = 1, . . . , 9. That is, the diagonal elements are 
J 2 , i = 1,2,4,5, while 



Ji • J2 — {J$ — Jf — ^2)/^ 

J2 • J5 = (^1 — J\ — Jt)/^ 



Ji ■ J4 — {Jj — J\ — Jf)/2, 

J4 • J5 = {j& — j\ — jI)/ 1 }- 



(5) 
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The two dot products that cannot be determined from the given lengths are u = Ji • J5 
and v — J2 - J4, which we regard as unknowns. These satisfy a linear equation obtained 
by squaring J 9 = -J 3 - J 6 , 

J 9 2 = J 3 2 + J 6 2 + 2{u + v + Ji • J 4 + J 2 • J 5 ). (6) 

Another equation connecting u and v is det G = 0, which holds since the four vectors 
lie in IR 3 and the 4-simplex defined by them is flat. This is a quartic equation in u and 
v, which by using to eliminate v can be converted into a quartic equation in u alone. 
We write this quartic as Q(u) = 0. We find the roots u of this quartic, solve for v by 
using ([6]), whereupon all components of the Gram matrix become known (there is one 
Gram matrix for each root). 

Ponzano and Regge (1968) discussed this procedure in somewhat different language, 
and apparently believed that all four roots would contribute to the asymptotics of the 9j- 
symbol. In fact, they do, if one wishes to work in the classically forbidden region and/or 
take into account tunnelling and exponentially small corrections in the neighborhood of 
internal caustic points (more about these below). But in the classically allowed region 
the asymptotics of the 9j-symbol are dominated by the contributions from "admissible" 
roots, namely, those roots that produce Gram matrices that can be realized as dot 
products of real vectors J«. Only these correspond to real geometrical figures of the 
type we have described. 

If a root u of Q(u) = is complex, then it produces a complex Gram matrix that 
cannot be realized with real vectors, and so u is inadmissible. But a real Gram matrix 
can be realized as the dot products of real vectors if and only if it is positive semidefinite, 
so even if u is real it will still be inadmissible if G has negative eigenvalues. 

We define the classically allowed region of the 9j'-symbol as the region in which 
Q{u) has at least one admissible root. In fact, in the classically allowed region Q(u) has 
four real roots of which two are generically admissible. We order the four real roots of 
Q(u) in the classically allowed region in ascending order and label them by k = 0, 1, 2, 3. 
It turns out that the two admissible roots are the middle two, k = 1,2, corresponding 
to the two terms of ([TJ with the same subscripts, k = 1,2. 

For a given admissible root, that is, a positive semidefinite Gram matrix, we wish 
to find the vectors Jj, i = 1,2,4,5. We arrange the four unknown vectors as the 
columns of a 3 x 4 matrix F, so that G = F T F. To find F given G, we diagonalize 
G, G = VKV T , where 1/ 6 0(4) and K is diagonal with nonnegative diagonal entries 
(the eigenvalues of G). At least one of these eigenvalues must be 0; we place it last, and 
write K = D T D where D is a real, 3x4 diagonal matrix. Then F = UDV T , where 
U is an arbitrary element of 0(3). This generates all possible sets of vectors whose dot 
products are realized in G; it amounts to using the singular value decomposition of F. If 
U = R e SO (3) then we generate a set of surfaces related by overall rotations; if U = — R 
we generate the time-reversed set. In this way a single Gram matrix, corresponding to a 
single admissible root of the quartic, produces a geometry and its time-reversed image. 
Altogether, the two admissible roots imply the four geometries in ([T]). 
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This method of finding F is discussed in the context of the 6j-symbol by Littlejohn 
and Yu (2009), where it is also applied in the classically forbidden region. There we find 
complex angular momentum vectors that satisfy the required algebraic relations. This 
carries over to the 9j-symbol in the classically forbidden region. In the literature on the 
6j-symbol it is common to state that a Euclidean group applies in the classically allowed 
region and a Lorentz group in the classically forbidden region; but for the 9j-symbol 
the groups are actually SO (3, R) and SO (3, C). 

5. The Classically Allowed Region and Configuration Space 

The classically allowed region is a subset of full dimensionality of the 9-dimensional 
parameter space of the 9j-symbol, itself a convex subset of M 9 defined by the triangle 
inequalities. To visualize this and other subsets of the parameter space it helps to fix 
seven of the j's to obtain a 2-dimensional slice. Figure H] illustrates such a slice for the 

case 



in which only j 3 and j$ are allowed to vary. The choice of j 3 and for this purpose is not 
arbitrary, since these two j's are quantum numbers for a pair of commuting operators 
on a space of 5-valent £77(2) intertwiners. They are like x and y for a wave function 
ijj{x,y). In this analogy, we think of (j3,je) -space as a "configuration space" for the 
9j-symbol and the 9j-symbol itself as a "wave function" We will mostly use 

the variables J 3 = j 3 + 1/2, J 6 = j 6 + 1/2 to describe this space. When thinking in 
classical terms, J 3 and Jq are continuous variables (not quantized). 

Figure 0] illustrates a convex region of the J 3 -J§ plane, bounded by straight lines 
and defined by the classical triangle inequalities, 



Properly speaking, configuration space is this convex region, not the whole plane. The 
unshaded area inside the convex region is the classically allowed region, surrounded 
by the shaded classically forbidden region. The caustic curve separates the classically 
allowed from the classically forbidden regions; it has kinks (discontinuities in slope) at 
points B, and is tangent to the boundary of the convex region at several points. Other 
features of this figure are explained below. 

Given a point (J 3 , J%) of the classically allowed region, the procedure described in 
Sees. |3] and H] produces a quartic polynomial Q(u) whose two middle roots k = 1,2 are 
admissible. These can be thought of as specifying a two-branched "root surface" that 
sits over the classically allowed region. The two middle roots coalesce as we approach 
the caustic curve, and become (inadmissible) complex conjugates as we move beyond. 




(7) 



maxfl Ji - J 2 |, I J e - J 9 |) < J 3 < min(Ji + J 2) J 6 + J9) 
max(| J 4 - J 5 |, I J 3 - J 9 |) < J 6 < min(J 4 + J 5 , J 3 + J 9 ). 



(8) 
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Thus, the two root surfaces can be thought of as being glued together on the caustic 
curve. 

Corresponding to each root there are two geometries modulo £0(3), related by 
time-reversal, so there is a two-fold "geometry surface" sitting above each root surface, or 
four geometry surfaces sitting above the classically allowed region. These four geometry 
surfaces are actually branches of the projection of an invariant 2-torus onto configuration 
space, and correspond to the four exponential terms in (pQ). This 2-torus sits in the phase 
space of the 9j-symbol, a 4-dimensional, compact symplectic manifold. 

This symplectic manifold is only one of several phase spaces that describe the 
classical mechanics of the 9j-symbol, but all the others have higher dimensionality so we 
call this one the "phase space of minimum dimensionality." It is one of the symplectic 
manifolds discovered by Kapovich and Millson (1996). Its analog in the case of the 
6j-symbol is a spherical phase space, which has been studied by Charles (2008) and 
by Littlejohn and Yu (2009). The phase space of minimum dimensionality is related to 
other phase spaces for the 9j-symbol by a combination of symplectic reduction (Marsden 
and Ratiu 1999) and the elimination of constraints. We have found it useful to employ 
all these spaces in our work on the 9j'-symbol. 



6. The Amplitude and Caustics 

The amplitudes of semiclassical approximations are notorious for the computational 
difficulties they cause. For example, several authors have resorted to computer algebra 
and/or numerical experimentation to check the amplitude determinant in the Ponzano- 
Regge formula. Actually, this amplitude (due originally to Wigner (1959)) is given by 
a single Poisson bracket between intermediate angular momenta (Aquilanti et al 2007, 
and, in more detail, Littlejohn and Yu 2009), which can be evaluated in a single line of 
algebra. More generally, semiclassical amplitudes are easily found in terms of matrices 
of Poisson brackets. 

In the case of the 9j-symbol we define 

Vijk = Ji ■ {Jj x J fe ), (9) 

which is six times the signed volume of the tetrahedron specified by edges i, j, k (it is 

the volume of the corresponding parallelepiped). Then the amplitudes Ax, A 2 in ([1]) are 
given by 

A = } (10) 

47r v /|detD| 

where 

D = ( Vu* V 215 \ (n) 

The subscripts 1,2 are omitted on A in ffTUj) because the same formula applies for both 
terms in (JTJ, but A\ ^ A 2 in general because the formula is evaluated on two different 
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geometries (associated with the two admissible roots). The quantitity detD is even 
under time-reversal, so the same amplitude applies to both a geometry and its time- 
reversed image. 

The volumes in matrix D are Poisson brackets of intermediate angular momenta 
in a recoupling scheme for the 9j-symbol, which are most easily evaluated in the phase 
space of minimum dimensionality. We omit details; suffice it to say for now that the 
derivation of the matrix ffTTl) in terms of Poisson brackets and thence the amplitude is 
extremely easy. 

We define the caustic set as the subset of the 9j-parameter space where det D = 0. 
Its intersection with the 2-dimensional slice seen in Fig. H] consists of the union of 
the caustic curve (the curve separating the classically allowed from the classically 
forbidden region) with the two points marked I. In addition, the caustic set includes the 
continuation of the caustic curves from points B into the classically forbidden region. 
The points I are "internal" caustics, that is, internal to the classically allowed region. 
While the caustic curve has codimension 1, the internal caustics have codimension 2. 

The quantity det D is nonzero away from the caustics. It turns out that the sign 
of det D distinguishes the two root surfaces, with det D > on root surface 1 and 
det D < on root surface 2. 

The caustics of the 6j-symbol occur at the flat configurations (flat tetrahedra), as 
appreciated by Ponzano and Regge (1968) and Schulten and Gordon (1975a,b). The 
caustics of the 9j-symbol, however, are not in general flat, that is, detD = does not 
imply that the configuration is flat. The flat configurations of the 9j-symbol, however, 
do lie on the caustic set. In a given J 3 - J 6 slice, there are precisely four flat configurations. 
In the example of Fig. HI these are marked B and /. The points B are flat configurations 
lying on the boundary of the classically allowed region (the caustic curve), while points 
I are internal flat configurations. As we vary the seven j's that are fixed in Fig. HJ 
the number of flat configurations on the boundary varies from 2 to 4; those not on the 
boundary are internal. 

In the usual manner of semiclassical approximations, (CQ) breaks down in a 
neighborhood of the caustic set (it diverges exactly at the caustic), and must be replaced 
by a diffraction function associated with a catastrophe (Berry 1976). In the case of the 
6j-symbol, the only catastrophe that occurs is the fold, yielding an Airy function as the 
semiclassical approximation, as noted by Ponzano and Regge (1968) and Schulten and 
Gordon (1975). This is the normal situation for systems of one degree of freedom. The 
9j-symbol, however, possesses two degrees of freedom, and other types of catastrophes 
occur. The fold catastrophe applies at most points along the caustic curve, where the 
9j-symbol is approximated by an Airy function; but at flat configurations there is an 
umbilic catastrophe, hyperbolic for those (B) falling on the boundary (caustic) curve and 
elliptic for the internal caustics (J). See Trinkhaus and Dreper (1977) for illustrations 
of the associated diffraction functions. The umbilic catastrophes are generic in systems 
of three degrees of freedom but occur in the 9j-symbol (with only two) because of time- 
reversal symmetry. However, only sections of the full three-dimensional umbilic wave 
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forms appear (Berry 1976). The cusp catastrophe, which can be expected in generic 
systems of two degrees of freedom, does not occur in the classically allowed region of 
the 9j-symbol. 

Caustics are associated with the coalescence of branches of the projection of a 
Lagrangian manifold in phase space onto configuration space. In the case of the 9j- 
symbol, the Lagrangian manifold is the invariant 2-torus mentioned in Sec. [5j Along 
the boundary of the classically allowed region, the two admissible roots coalesce, which 
means that the four geometries merge into two. At most points on the boundary 
curve, the two remaining geometries are not equal, but are related by time-reversal. 
At such points we have a fold catastrophe, and the 9j-symbol is approximated by an 
Airy function (modulated by a cosine term). At points B, however, the two geometries 
related by time-reversal merge into a single flat configuration, producing the hyperbolic 
umbilic catastrophe. 

At internal caustic points / the geometry and its time reversed image for one of 
the two admissible roots coalesce to produce a flat configuration. The two geometries of 
the other root surface, however, do not coalesce. Thus at internal caustics I there are 
three geometries. Only the flat configuration associated with one of the roots produces 
the elliptic umbilic catastrophe; thus, only one of the two terms in ([1]) is replaced by 
the elliptic umbilic diffraction function, while the other remains as shown in ([1]). The 
9j symbol is a linear combination of these two terms, but the elliptic umbilic diffraction 
function dominates when the scaling factor k is large. 

The caustics have a certain size, that is, a distance around the caustic set over 
which diffraction functions must be used instead of ([IT) . This distance Aj scales as k 1 ^ 3 
for all three catastrophe types (fold and elliptic and hyperbolic umbilic) discussed here. 

In the neighborhood of fold catastrophes the wave function scales as k~ 17 ^ 6 , that 
is, k l l & higher than the k~ 3 of the two terms in ([1]). In the neighborhood of umbilic 
catastrophes the scaling is k~ 8 ^ 3 , that is, with another factor of fc 1//6 . For large values 
of k the 9j-symbol is largest near the points /, B. 

Linear combinations with different scaling behaviors have been observed by Barrett 
and Steele (2003) and by Freidel and Louapre (2003) in their studies of the lOj-symbol. 
It seems that the 9j-symbol is the simplest spin network in which this phenomenon 
occurs. 

7. The Phase 

The phases S\ and S 2 in ([Tj) each have the form 



where 6i is the angle between normals of adjacent faces of the geometrical figure. This 
of course is similar to the Ponzano-Regge formula, but the 6j-tetrahedron is convex 
and all dihedral angles can be taken in the interval [0,7r]. The dihedral angles for the 
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(12) 



i=l 
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9j-symbol, on the other hand, must be allowed to lie in a full 2ir interval, as explained 
momentarily. The subscripts 1,2 are omitted on S in (fT2l because the same formula 
applies to both terms in (TjQ). The formula must be evaluated, however, on two different 
geometries, so S\ and S 2 are not equal. In addition, the angles 9i lie in different intervals 
for the two geometries. 

Each edge i of the geometrial figure is adjacent to two faces, for example, edge 4 in 
Fig. [2]is adjacent to faces l'A'T and 456. One face adjacent to edge i contains vector Jj, 
and the other J[. Let the two normals of these two faces, according to the conventions 
given above, be n and n'. Then we define 9i as the angle such that 

R(j,ei)n = n', (13) 

where j is the unit vector along J, specifying the axis of a rotation R by angle 9i using 
the right-hand rule. In the Ponzano-Regge formula one can compute the dihedral angle 
from its cosine, but for the 9j one must also use the sine of the angle. That is, ( fTBl) is 
equivalent to 

n' = cos 9i n + sin 9j j x n. (14) 

This determines 9{ to within an additive integer multiple of 2n. We add the further 
requirement that for the geometries associated with the first root (the cosine term in 
p])), — 7r < 9i < +7i, while for the second root (the sine term in (JTJ)), < 9i < 2-n. These 
ranges for the angle 9{ are chosen because they give a continuous branch for the angle 
over the two root surfaces. It turns out that 9i never crosses ±7r on the surface for root 
1, and it never crosses or 2n on the surface for root 2. 

The rules given in Sees. |3] and H] for converting vectors into surfaces with oriented 
edges and triangles are an essential part of the definition of the dihedral angles 9{. It is 
of interest to see how the angles change when a set of vectors or the associated geometry 
is subjected to some symmetry. 

Under time-reversal, the orientation of all triangles reverses, that is, the normal 
vectors stay the same but the vectors defining the edges are inverted. This means that 
the angles 9i go into —9{ on root surface 1, while they go into 27r — 9{ on root surface 2 
(both changes guarantee that the angles remain within their respective ranges). Thus 
S goes into —S on root surface 1 and 

S -»• -3 + 2nv + 9?r (15) 

on root surface 2, where v is the integer, 
9 

" = X> (16) 

i=i 

These guarantee that cos 3\ and sin S2 are invariant under time-reversal. Since the 
same applies to the amplitudes A\ and A 2 , one can choose either a geometry or its 
time-reversed image, for each root, when evaluating (CQ). 

This completes the definition and geometrical interpretation of all the notation used 
in ©. 
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8. Symmetries of the 9j-symbol 

The formula (JTJ) transforms correctly under the symmetries of the 9j-symbol 
(Varshalovich et al 1981, Sec. 10.4), which state that the 9j-symbol suffers a phase 
change of (— l)" under odd permutations of rows or columns or under transposition. 
Consider, for example, the swapping of the first two columns, and let P be the 
permutation of indices, so that PI = 2, P2 = 1, P3 = 3, etc. This maps an old 
set of nine j's into a new set, and old quartic Q(u) = into a new one, etc. We find 
that the u root of the old quartic becomes the v root of the new one, which amounts to 
saying that the root 1 surface of the old geometry is mapped into the root 2 surface of the 
new one, and vice versa. Also, the orientations of the three unprimed triangles reverse, 
but not those of the primed ones, causing all nine dihedral angles to be incremented or 
decremented by tt (depending on the range). If we let 9j be the original angles and 9{ 
the new ones, then when 9i is on root surface 1 we find 9pi = 9i + tt, which means that 
the new angle is in the right range since it is on root surface 2. Similarly, when 9\ is on 
root surface 2 then 9pi = 9{ — tt, which is in the right range since 9pi is on root surface 
1. As a result, when the original geometry is on root surface 1, we have 



so that sin ^2 = ( — 1)" c °s Si, while if the original geometry is on root surface 2, we have 



so that cos S\ = (— iy sin 5*2. The sine and cosine terms in (JTJ swap under column swap, 
and the result acquires an overall phase of (—l) u , as required. The specified ranges on 
the dihedral angles on the two root surfaces are necessary for this to work out. 

9. Comments and Conclusions 

It is easy to derive the expression (jl2l) by the method of Roberts (1999), which involves 
rotating faces by an angle of tt about their normals, and edges by an angle of tt about a 
normal to them. The phase ( jl~2l) (times 2) is then an action integral along one Lagrangian 
manifold and back along another (the analogs of the A- and 5-manifolds of Aquilanti 
et al (2007)). Similar expressions apply to any spin network of any complexity. But the 
contours chosen for the integration are not unique, in that one can add any multiples of 
quantized loops on the two manifolds. These modify both the actions and the Maslov 
indices, and amount to changing the choice of branch for the angles 0$, that is, adding an 
integer multiple of 2tt to these angles. This does not leave the trigonometric functions in 
(CQ) invariant because the angles are multiplied by the J{, which may be half-integers. The 
result is that the phase of the approximation to the 9j-symbol depends on the contours. 
A more serious worry is that the contours, that is, the branches for the 9i, may change 




(17) 



i=i 



i=i 




(18) 
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as we move around in the parameter space of the 9j-symbol. This would amount to 
crossing a branch cut for the angles 9i (and there are different branch cuts for different 
angles). In addition, as we move around in parameter space we can make any two 
adjacent faces rotate relative to one another around their common edge as many times 
as we want. Although the phases in question are "only" powers of —1, straightening 
out this issue was by far the hardest part of this work. In the end we realized that the 
ranges [— n, +7r) on root surface 1 and [0, 2tt) on root surface 2 guarantee that there are 
no branch cuts and hence no discontinuities. The ranges specified for the angles d{ give 
us in effect a global, smooth definition of contours for carrying out action integrals. 

We present several numerical comparisons of (JTJ) with the exact 9j-symbol. In 
Fig. [5] the approximation ([!]) (smooth curve) may be compared to the exact 9j-symbol 
(sticks) as a function of ja for fixed values of the other j's. The range chosen lies inside 
the classically allowed region, far from a caustic. Fig. [6] shows the comparison in a 
range that crosses a fold catastrophe, and Fig. [7] shows the comparison in an interval 
that passes near a hyperbolic umbilic catastrophe (the upper point I in Fig. HJ). The 
approximation (TjQ) is too large near the point /. 

Varshalovich et al (1981) present an asymptotic approximation for the 9j-symbol 
without citation (their Eq. (10.7.1)). We believe their formula must be an asymptotic 
expression for the 9j-symbol in a different sense than we have defined it, or else it is 
incorrect. 

The two terms in (pP) have different trigonometric functions (sine and cosine) 
because there is a relative Maslov index of 2 between the two root surfaces. The 
relative Maslov index between a geometry and its time-reversed image is 0, a somewhat 
surprising result because in mechanical systems and in the 6j-symbol the Maslov index 
between a branch or geometry and its time-reversed image is 1. 

When an interior caustic occurs on a root surface, the two geometries that sit 
above it form a double cover, in the manner of the Riemann sheet for the square root 
function. The internal caustic point I is a branch point for the cover. Geometries 
transform continuously into their time-reversed images as we go around the point /, 
without crossing a caustic. 

Several studies of the asymptotics of spin networks have started with an integral 
representation of the network, to which the stationary phase approximation is applied. 
Roberts (1999) represented the 6j-symbol as a scalar product in a certain Hilbert 
space, which was put into the coherent state representation, whereupon the integral 
was evaluated by the stationary phase approximation. Coherent states have played a 
prominent role in many recent semiclassical studies. Our approach has been to work as 
much as possible in a representation-independent manner. For example, the stationary 
phase points are seen as intersections of Lagrangian manifolds. Some of the basics of 
this approach were presented in Aquilanti et al (2007). We have not specifically used 
the coherent state or any other representation. 

Some aspects of this calculation carry through in an obvious way to higher spin 
networks, while for others nontrivial generalizations seem to be required. But we believe 
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that an understanding of the 9j results are necessary for a full understanding of the 
asymptotics of higher spin networks. 

We will report in more detail on the derivation of (CEJ) in a later publication. 
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Figure 1. The spin network for the 9j-symbol. 



Figure 2. The six triangles defined by (0| form the "central region" of the figure, 
with three triangles duplicated and attached to the edges of the central region. The 
notation 1, 2', etc refers to Ji, 3' 2 , etc. 



Figure 3. The central region of the time-reversed surface. Notation —1, —2' etc refers 
to — Ji, — 3' 2 , etc. 



Figure 4. The convex region of the J^-Jq plane is the configuration space of the 
9j-symbol. The shaded area is the classically forbidden region, and the unshaded, the 
classically allowed. Points I are internal caustics, two of the flat configurations; points 
B are the other two flat configurations, lying on the boundary curve. 



Figure 5. Comparison of exact 9j-symbol (vertical sticks) with approximation |T]), 
away from a caustic. Values used are those in ([7|), with j$ = 50. 



Figure 6. Like Fig. [H but an interval that spans a fold catastrophe (with j$ = 60). 
The approximation fl} is discontinued at the caustic, the exact values are continued 
into the classically forbidden region. 



Figure 7. Like Fig. [SI but passing near an ellipitc umbilic catastrophe (with j 6 = 79). 



